Towards a fully consistent parameterization of modified gravity 
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There is a distinct possibility that current and future cosmological data can be used to constrain 
Einstein's theory of gravity on the very largest scales. To be able to do this in a model-independent 
way, it makes sense to work with a general parameterization of modified gravity. Such an approach 
would be analogous to the Parameterized Post-Newtonian (PPN) approach which is used on the scale 
of the Solar System. A few such parameterizations have been proposed and preliminary constraints 
have been obtained. We show that the majority of such parameterizations are only exactly applicable 
in the quasistatic regime. On larger scales they fail to encapsulate the full behaviour of typical models 
currently under consideration. We suggest that it may be possible to capture the additions to the 
'Parameterized Post-Friedmann' (PPF) formalism by treating them akin to fluid perturbations. 



I. INTRODUCTION 

It is possible that we live in a Universe in which more 
than 96% of the energy and matter density is in the form 
of an exotic dark substance. The conventional view is 
that roughly a quarter of this obscure substance is in the 
form of dark matter and the remainder is in the form 
of dark energy. Theories abound that propose explana- 
tions for dark matter and dark energy and there is an 
active programme of research attempting to understand 
and measure them. 

It may also be possible that our understanding of grav- 
ity is lacking, and that Einstein's theory of General Rela- 
tivity (and more specifically, the Einstein field equations) 
are not entirely applicable on cosmological scales. The 
past decade has seen unprecedented growth, from a hand- 
ful to a veritable menagerie of possible modifications to 
ravity that may be perceived as a fictitious dark sector 

The proliferation of theories of modified gravity 
couldn't have come at a better time. Observational cos- 
mology has entered what some have called an era of 'pre- 
cision cosmology'. Hubristic as such a point of view 
might be, it is certainly true that cosmology is being 
inundated by data, from measurements of the Cosmic 
Microwave Background (CMB) [3, Qj galaxy surveys J4l]i 
weak lensing surveys [5| and measurements of distance 
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and luminosity at high redshift with supernovae la j(|. 
With such data in hand it is possible to test cosmolog- 
ical models and constrain their parameters with some 
precision. With the forthcoming experiments currently 
on the drawing board Q-Q' great things are expected. 
In particular, there is a hope that it may be possible to 
distinguish between the two paradigms: the dark sector 
versus modified gravity. 

The situation in cosmology is reminiscent of that in 
General Relativity in the late 1960's and early 1970's. 
Then, Einstein's theory was undergoing a golden age with 
discoveries in radio and X-ray astronomy, as well as pre- 
cision measurements in the Solar System and beyond, 
making it increasingly relevant. As a result, a plethora 
of alternative theories of gravity were proposed which 
could all in principle be tested (and ruled out) by ob- 
servations |10l - [l3 |. Out of this situation a phenomeno- 
logical model of modified gravity emerged, the Param- 
eterized Post- Newtonian (PPN) approximation J13rtl6| , 
which could be used as a bridge between theory and ob- 
servations. In other words, from observations it is possi- 
ble to find the constraints on the parameters in the PPN 
approximation. The constraints are model-independent. 
From any given theory it is then possible to calculate the 
corresponding PPN parameters and find if they conform 
to observations. The PPN approximation is sufficiently 
general that it can encompass almost all modified theo- 
ries of gravity that were then proposed. 

Clearly something like the PPN approach is desirable 
in cosmology. Given the rapid increase in the number of 
modified theories of gravity, it would make sense to con- 
struct a parameterization that could serve as a bridge be- 
tween theory and observation. Observers could express 



their constraints in terms of a set of convenient parame- 
ters; theorists could then make predictions for these pa- 
rameters and check if their theories are observationally 
viable. Instead of performing many constraint analyses 
on individual theories, one could run just a single con- 
straint analysis on the parameterized framework. With 
a dictionary of translations between theories and the pa- 
rameterization in hand, these general constraints could 
be immediately applied to any particular theory. 

Another key advantage of a parameterized approach is 
that it allows one to explore regions of theory space for 
which the underlying action is not known. For example, 
in WIBI we will see how the Lagrangian f{R) is related to 
our framework 'parameters' (which are really functions, 
not single numbers - see below). Cosmological data will 
exclude certain regions of parameter space. If a new form 
of f{R) is proposed in the future, it should be a quick 
operation to see whether it falls into the excluded region 
- even though that particular Lagrangian was not known 
at the time that the constraint analysis was performed. 

In this paper we will discuss the requirements of how 
to parametrize modifications to gravity on cosmological 
scales. It develops the principles first put forward in [17[ 
and explores how they may be applied more generally. 
The layout of this paper is as follows. In S|TT] we discuss 
the idea behind the PPN approach and show that it can't 
be imported wholesale into cosmology. We then briefly 
look at attempts at parametrizing gravity and point out 
their limitations. In §1111 we discuss a possible formal- 
ism in detail - we co-opt the name Parameterized Post- 
Friedmann approach from 18| - and argue that it may 
be sufficiently general to encompass a broad class of the- 
ories. In §IVI we construct the hierarchy of equations 
that should be satisfied in the case in which there are no 
extra fields contributing to the modifications to gravity. 
In fjV] we discuss the more general case with extra fields 
and how this affects the relations between the different 
coefficients. In WII we focus on four modified theories of 
gravity and analyse how they fit into the formalism that 
we are proposing. Finally, in Will we summarize the 
state of play of the parameterization we are proposing. 



<&) are gravitational potentials and 7ppn is a PPN param- 
eter, equivalent to one of the older Eddington-Robertson- 
Schiff parameters. 

There are a few properties which are of note in this ex- 
pression. First of all, the parameterization is constructed 
around a solution of the Einstein field equations, the 
Schwarzschild solution (with 7ppn = 1) - the field equa- 
tions do not come into play. Second (and this isn't obvi- 
ous from the expressions above), the parameter 7ppij only 
depends on parameters in the theory and not on inte- 
gration constants or 'environmental' parameters such as 
the central mass. This means that, given a theory, it 
is possible to predict 7ppn solely in terms of fundamen- 
tal parameters of the theory (i.e. the parameters in the 
action). Finally, we see that the mismatch between the 
gravitational potentials can be expressed as [iJl 



$ - * = Cppn* 



(2) 



with Cppn = (7ppn — 1)/7ppn often called the gravitational 
slip. General Relativity is recovered when Cppu — 0. 

The idea of applying an equation of the form of eqn. ([2]) 
to cosmology emerged from the work of Bertschinger in 
[iSi . Bertschinger showed that on large scales it was pos- 
sible to calculate the evolution of ^ and $ using only in- 
formation about the background evolution and assuming 
a closure relationship between the two potentials. The 
simplest assumption is a closure relation of the form of 
eqn. ([2]), but in no way was it implied that this would be 
a realistic relationship that would be valid in the general 
space of theories of modified gravity. 

Nevertheless, over the last few years the simplified 
equation for gravitational slip has been adopted as a gen- 
eral parameterization which should be valid in cosmology 
[20:- 22]. It has been shown to be valid in a few cases, in 
the quasistatic regime (i.e. on small scales), and explicit 
expressions have been found for C in terms of fundamental 
parameters of those theories (some examples are collected 
in 1]). Such a parameterization has been extended to in- 
clude another parameter, a modified Newton's constant 
Gcff, which may differ from Gq. The method is then to 
use eqn.([2|) and a modified Poisson equation, 



II. THE CONVENTIONAL APPROACH 



-fc2$ = 47rGcff a^ V PiAi 



(3) 



The plan is to construct a parameterization that might 
mimic the PPN approach on cosmological scales. It 
is therefore useful to look very briefly at the PPN ap- 
proach, which proceeds as a perturbative expansion in 
v/c (though we set c=l in what follows). Consider the 
modified, linearized Schwarzschild solution: 



ds^ 



(1 - 2*) dt^ + (1 - 2$) dr^ + r^dn^ 



GqM 

r 
7ppn* 



(1) 



where Go is Newton's constant, M is the central mass, fi 
is the two-dimensional volume element of a sphere, (^, 



(where pi is the energy density of fluid i and A^ is the co- 
moving energy density) to modify the evolution equations 
for cosmological perturbations. A modified Einstein- 
Boltzman solver is then used to calculate cosmological 
observables. The two parameters (GcS, C) have been 
adopted more generally and have been used to find pre- 
liminary constraints on modified gravity theories by a 
number of groups |2l| - [24j . 

Clearly such an approach to parametrizing modified 
gravity has some significant differences with the PPN 
approach. For a start, modifications are applied to the 
field equations and not to specific solutions of Einstein's 
field equations. This is understandable - the solutions 



of interest in cosmology are not only inhomogeneous 
but time-varying, unlike the incredible simplicity of the 
Schwarzschild solution that arises due to Birkhoff 's the- 
orem. Also, unlike in PPN, the parameters at play - ( 
and Goff - will not only depend on fundamental param- 
eters of the theory but also on the time evolution of the 
cosmological background. 

Ideally, any time-dependence in the parameterization 
will be simply related to background cosmological quanti- 
ties (like the scale factor, energy densities or any auxiliary 
fields that are part of the modifications) and not depen- 
dent upon the time evolution of $, ^ or any other per- 
turbation variables. For such a requirement to be possi- 
ble it is essential that any parameterization is sufficiently 
general to encompass a broad range of theories. As we 
will show in this paper, parameterizations using eqns.Q 
and ^ are simply not general enough to capture the full 
range of behaviour of modified theories of gravity. It has 
been argued that such a parameterization can be used as 
a diagnostic; that is, for example, a non-zero measure- 
ment of C might indicate modifications of gravity |25| . 
This may be true, but such a measurement cannot then 
be used to go further and constrain specific theories. It 
would be more useful to build a fully consistent parame- 
terization which can be used as a diagnostic and can be 
linked to theoretical proposals. The purpose of this paper 
is to take the first steps towards such a parameterization. 



III. THE FORMALISM 

When considering modified gravity theories it can be 
helpful to cleanly separate the non-standard parts from 
the familiar terms that arise in General Relativity (hence- 
forth GR). We can always write the modifications as a 
an additional tensor appearing in the Einstein field equa- 
tions, i.e. 



Gfj,iy — SirGoa T^^ + a L/^i, 



(4) 



The diagonal components of the tensor U^u are equiva- 
lent to an effective dark fluid with energy density X and 
isotropic pressure Y (where the constants have been ab- 
sorbed). The zeroth-order Einstein equations are then: 

Ef = ZT-L^ + iK^^TiGaa^^Pi + a^X (5) 

i 

Er = -{2H + H^ + K)^8TTGoa^J2Pi + "-^^ (6) 

i 

where H — H/a is the conformal Hubble parameter and 
K is the curvature. We will use Ep and Er as defined 
above throughout this paper. For future use we define 
E = Ep + Epj. The summations in the above expressions 
are taken over all conventional fluids and dark matter, 
and dots denote derivatives with respect to conformal 
time. In this paper we will largely adhere to the defini- 
tions and conventions used in [l3| . We also note that the 



Bianchi identity V^G(^ — implies the relation: 
Ep+niEF + SER) = 



(7) 



Assuming that the conservation law V ^Tj^ — holds sep- 
arately for ordinary matter and the effective dark fluid, X 
and Y must be related by the equation X-t-3'H(Ar + y) = 
0. 

Gontinuing in this vein, our goal is to write the linearly 
perturbed Einstein equations as: 



6G^y = ^ttGqo^ST^i, + a^5U^y 



(8) 



In general the tensor SU^j^n will contain both metric per- 
turbations and extra degrees of freedom (hereafter d.o.f) 
introduced by a theory of modifled gravity. We can sep- 
arate 5Uf^i, into three parts: i) a part containing only 
metric perturbations, ii) a part containing perturbations 
to the extra d.o.f., iii) a part mixing the extra d.o.f. and 
perturbations to the ordinary matter components: 



5U^, 



(9) 



The argument variables in this expression will be intro- 
duced shortly. We have written the Einstein field equa- 
tions such that Tj^ contains only standard, uncoupled 
matter terms and hence obeys the usual (perturbed) con- 
servation equations, 5{y ^T^}) = 0. As a result U^ must 
obey its own independent conservation equations, so that 
at linear order we have 5{V ^U^^) = 0. We will use the fol- 
lowing notation to denote components of SUj^ from here 
onwards: 



Up = a'dUl D,jU^ - a^(U; - ^U^S}) 



(10) 



where Dij — V^Vj 



l/3gy V^ projects out the longitu- 



dinal, traceless part of 5Uj^ and qij is a maximally sym- 
metric metric of constant curvature K. The definition 
of Us. In the case of unmodified background equations 
perturbed conservation equation for Uj^ gives us the fol- 
lowing two constraint equations at the linearized level 

Ua + nUA - V^C/e + ^a^{X + Y){$ + 2^^e) 

+HUp = 0(11) 
Ue + 2nUe-^Up-^V^Us + a\X + Y)E = 0(12) 

where the metric fiuctuations /3, e and S are defined in 
eqn. (fT3| . 

In this paper we will initially present general forms for 
the construction of metric-only ^C/^i, that satisfy equa- 
tions ([TT|) and ([T2|). then impose the restriction that 
the field equations can contain at most second-order 
derivatives. Gravitational theories containing derivatives 
greater than second-order are generally disfavoured as 



they typically result in instabilities or the presence of 
ghost solutions [ll, [2^[23|- However, we note that some 
special cases of higher-order theories are acceptable e.g. 
f{R) theories [28| (see WIBI) . Hence we start with the 
general case in order to indicate how our results may be 
extended to higher-derivative theories |29l |. 

The requirement of second-order field equations means 
that Ua and Uq can only contain first-order derivatives 
with respect to conformal time, as can be seen from 
egns. pTj ) and (|12p . The specific implications this has 
depends on which of the tensors in eqn.® are present. 
In mV\ we will explore the structure of theories with only 
metric perturbations, whilst theories with extra degrees 
of freedom will be presented in jVl In Appendix |B] we 
display formulae for generating constraint equations in 
an arbitrary-order theory of gravity with no additional 
degrees of freedom. 

We write the perturbed line element (for scalar pertur- 
bations only) as: 

ds^ = - a'^{l-2E)dt'^ -2a^{V,e)dtdx 



+ a^ 



'"--3^ 



q^J 



D,. 



dx* dx' 



(13) 



It will prove useful to define the gauge- variant combina- 
tion: 

V^v + 2e (14) 

We also define the gauge-invariant potentials: 



$ 



* 



(/3 



\v-\'HV 



(15) 

$ and 'I' are equivalent to the Bardeen potentials — ^h 
and ^A respectively. Note that ^ contains a second- 
order time-derivative. In the first four sections of this 
paper we will frequently use a linear combination of these 
variables that remains first-order in time derivatives of 
perturbations, due to a cancellation between the V terms: 



(17) 



From f|V] onwards we specialise to the conformal Newto- 
nian gauge, and hence revert to the familiar potentials 
<I> and \&. We introduce a shorthand notation for the 
components of SG^v in exact analogy to that introduced 
for SU^v: i.e. -Ba = —a^SG^ etc. Hereafter the left-hand 
sides of the perturbed Einstein equations will be denoted 
by: 



Ea = 2{V^ + 3K)<^ 
Eq = 2 fcf + -EV 



enkv - -HEV 



Ep = &k—- 
dr 



WHkT - 2{V^ + 3X)($ - *) 



Ey 



-3£;* + - {Eb,-2HEr]V 



In terms of these variables the perturbed Einstein equa- 
tions are Il7t: 



Ea = 8^Ga2^p,(5, + C/a 

i 

Ee = 8^Ga2^(p, + P0^» + C^e 

i 

Ep = 247rGa2 ^ p.Hj + C/p 

i 

Es = 8^Ga2^(p, +P,;)S^ + C/E 



(19) 
(20) 

(21) 
(22) 



For simplicity we will hereafter consider only the case of 
a universe with zero spatial curvature, K — 0. 



IV. THE GENERAL PARAMETERIZATION - 
NO EXTRA FIELDS 

A. General case - unmodified baclcground 

Let us begin with the simplest case by applying two 
restrictions: i) We consider the case of modifications to 
gravity that appear only at the perturbative level, that 
is, they maintain the background equations of GR for 
a Friedmann-Robertson- Walker metric; ii) there are no 
new d.o.f. present in the theory, so SUfj^i, contains only 
metric perturbations. We will relax restriction i) in ijIVDI 
and restriction ii) in ^V) We will see shortly that the 
treatment presented in this subsection is also applicable 
to ACDM, because the X + Y terms in eans.([TT|) and 
(IT^ vanish for a cosmological constant. The agreement 
between an exact ACDM background and current data 
means that theories obeying restrictions i) and ii) are 
of particular interest, even though they correspond to a 
limited region of theory space. 

The requirement of gauge form-invariance places 
strong restrictions on the forms that 6U^ij can take jl7| . 
We will postpone a detailed discussion of these restric- 
tions until ijIVDl where they will be a useful tool in guid- 
ing us to allowed combinations of metric perturbations. 
In this subsection it suffices to point out that the stan- 
dard Einstein field equations of GR are of course already 
gauge form-invariant; so any additive modification like 
SU^i, must be independently gauge form-invariant in or- 
der to preserve the invariance of the whole expression. 
This property is a direct consequence of the fact that we 
have not yet modified the background equations. In this 
case the only objects that can be present in the tensor 
6U^„ are the gauge- invariant metric potentials $ and F. 

So, we can construct the tensor SU^i, from series of 
all the possible derivatives of $ and F. This structure 
should be general enough to encompass any metric the- 
ories, where the action is constructed purely from curva- 
ture invariants, e.g. /(i?) gravity, Gauss-Bonnet gravity 
|30| and Lovelock gravity [3l| . If we wish to parameterize 
(18) only second-order theories then we will need to truncate 



these series at A^ = 2, as discussed in section IIIII The 
components of t/^^ are given by: 



^-2 



n=0 

N-1 

n=0 
N-1 



Tl = 



(23) 
(24) 
(25) 
(26) 



The coefficients An-Kn are functions of the scale factor 
a, wavenumber k and background quantities such as pi - 
for the sake of clarity we will suppress these dependencies 
throughout. The factors of k ensure that the coefficient 
functions are dimensionless. 

Let us take a moment to explain the upper limits on 
the summations in eans. (P5) )- (l^^ . $ and F are first-order 
in time derivatives (see eqns. p3|) and ([TT])). Ua is differ- 
entiated in egn. dTTI) . so truncating the series in eqn. ((23)) 
at $(^~2) gives field equations containing time deriva- 
tives of order N. Uq is treated analogously to Ua ■ As Up 
and Us are not differentiated in the components of the 
Bianchi identity, the series in eans. (p5) ) and (|26)) are al- 
lowed to extend one order higher than those in eqns.([23]) 
and (Pll) . 

We substitute our forms for Ua, Uq, Up and t/s hi^to 
the components of the Bianchi identity ([TT|) and p2)) . $ 
and r are non-dynamical fields and so will not evolve 
in the absence of source terms. Yet when we perform 
the substitution, the Bianchi identity appears to give us 
two evolution equations for $ and F. The only way this 
can be avoided is if the coefficients of each term 'l'^"^ 
and F^"^ vanish individually, which provides us with con- 
straint equations on the functions An-Kn (this procedure 
will be clarified with an example shortly). Each compo- 
nent of the Bianchi identity results in N constraint equa- 
tions from each of the $'"' terms and F'"^ terms, and 
eans. (P5) )- (PB| contain 8N-4 coefficient functions in total. 
Hence we have 4N — 4 free functions with which to de- 
scribe the theory. 



B. Second-order case - unmodified bacl^ground 

In Appendix |B] we give formulae for generating the 
constraint equations of an arbitrary-order theory with 
unmodified background equations. We will now explic- 
itly present the second-order case, which corresponds to 
setting N=2 in eans. (P5)) - (l^ . In a general case this will 
give us four free functions. However, if the background 
equations are unaltered then we must set Fq = Iq = be- 
cause F contains a second-order conformal time deriva- 
tive of the scale factor. One might consider cancelling 



this a term by adding a term proportional to EV, but 
this would break the gauge-invariance of the perturbed 
Einstein equations. We will see later that modification of 
the background equations allows us to add an EV term 
without violating gauge-invariance, which in turns means 
that F can be present in Ua and Uq. 

Using eqns. ([m) and (IT^ we find that setting Fq = Iq = 
forces Ji — Ki — also. Then, for the second-order 
case, the remaining terms in 5Uf^y are: 



Ua 
Ue 
Up 






Jofc^f 



= Do<i> + ^<i> + KaT 



(27) 



The constraint equations are given in Table HI indicat- 
ing the terms and Bianchi identity from which they arise 
{Bl ^ egn. pTj) . B2 ^ eqn. (|12p V These expressions can 
be generated using the formulae in Appendix |B] We can 
see immediately that the F terms in Up and C/s vanish, 

leaving SUfj,^, expressed entirely in terms of $ and $. We 
have two free functions remaining, which we will choose 

to be Dq and Di. Eliminating Ci from the two <t con- 
straints gives (where Hk = H/k): 



jD, = -^iAo + 3HkBo) 



(28) 



The combination on the right-hand side appears when we 
form the (Fourier-space) Poisson equation from eqns. p^ 
and (|20p . where it acts to modify the value of Newton's 
gravitational constant: 



-fc2$ 



4:TT- 



( 

1-3 



^^E 



Pi^l 



(29) 



where Ai =,5, -I- 3^(1 
perturbation and 



1 



9 



{Ao + 3HkBo) 



-Wi)6i is a gauge-invariant matter 



(30) 



The sum in eqn. (j29p is over all known fiuids and dark 
matter, and Go denotes the canonical value of Newton's 
constant. From here on we will replace Di/k in eqn. (|27|) 
by g/l-L to remind us of the connection between the mod- 
ifications to the slip relation and the Poisson equation. 
We will also replace Dq by C, to distinguish it from the 
other coefficient functions, which can all be expressed in 
terms of g and C, using the constraint equations. We con- 
tinue to suppress the arguments of g and C- 

The effective gravitational constant appearing in the 
Poisson equation is Geff = Go/(l — g). The traceless 
space-space component of the Einstein equations be- 
comes: 



$_>jf 



BTrGo^(pj 



Pi)^: + C* - 



^4 



(31) 





Origin 


Constraint equation 


1 


[Bl] $ 


io + nAo + kBo + nCo = 


2 


[Bl] $ 


^0 + WfcCi = 


3 


[Bl] f 


Jo = 


4 


[B2] $ 


Bo + 2HBo - ifcCo + 2fcDo = 


5 


[B2] $ 


Bo - iCi + iDi = 


6 


[B2] f 


27^0 - Jo = 



TABLE I: Table ol the constraint equations for the second- 
order metric theory specified in ijlVBI These can be generated 
using the formulae in Appendix [B] 



In parameterizations equivalent to {Q, rjsUp) the degener- 
acy is Q{1 + l/jysiip), so for lensing applications it makes 
sense to define new parameters along and perpendicular 
to the degeneracy direction |32434j . In the {g, () parame- 
terization a degeneracy remains. The dominant contribu- 
tions to lensing signals come from quasistatic scales, on 
which time derivatives of perturbations can be neglected 
(see later for a fuller discussion). The degeneracy is then: 



-fc2($ + f) = AnGoa^ ^ PiS, 



(2-C) 



(34) 



The anisotropic stress perturbation E; is automatically 
gauge- invariant, but negligible for standard fluids at late 
times. The above expression echoes its PPN equivalent, 
eqn.(l2]); but note that, as discussed in S|lTl (^ is a func- 
tion of background quantities (which potentially intro- 
duce time- and scale-dependence), whereas Cppn depended 
only upon fundamental parameters of a gravitational the- 
ory. 

Other authors have made numerous different choices 
for the two free functions of a second-order theory; a 
useful summary of some of these is provided by [22| . A 
common choice is to introduce a function Q = Gcs/Go, 
related to our g hy Q — {1 — g)~^ [21| (though dif- 
ferent notation is in no short supply) . The relation- 
ship between the two potentials is often parameterized 
as $ = T?siip(«i k) ^ in the spirit of the PPN parameter 
IPPN [13']. It might be felt that by introducing yet an- 
other parameterization of PPF we are adding to this dis- 
array. However, in the next subsection we will argue that 
a two- function slip relation such as eqn. (1311) is needed to 
avoid implicitly introducing higher-order derivatives into 
a purely metric theory. 

Writing the relationship between the two gauge- 
invariant potentials as $ = •qs\ip{a,k)'^ implies that the 
spatial off-diagonal component of the Einstein field equa- 
tions is: 



$-* = ( 1 I $ 



(32) 



Comparing the above equation with eqn. pip implies: 



'7siip 



= i-c- 



g dln$ 
n dT 



(33) 



Now 77siip has an environmental dependence, which is 
problematic. We would require detailed knowledge of the 
environment in which we wish to test a theory a priori, 
and the PPF functions would need to be recalculated for 
numerous different situations. Unless <& = 0, the param- 
eterizations in egns. ipTj ) and ([5^ do not have a simple 
equivalence. 

A degeneracy arises between g and C when comparing 
to data from weak gravitational lensing, which probes the 
combination $ -f ^If in the conformal Newtonian gauge. 



It seems that neither of the two parameterizations pre- 
sented so far are optimal for weak lensing constraints. 



Why neglecting g in the slip relation implies a 
higher-derivative theory 



We have seen in the previous section that in a metric- 
based second-order theory of modified gravity the most 
general form of the gravitational slip should be expressed 
in terms of the gauge-invariant potential $ and its first 
derivative with respect to conformal time. Two free func- 
tions C and g were used as the coefficients of these terms 
respectively, where g resulted in a modification to New- 
ton's gravitational constant in the Poisson equation. Us- 
ing a single function to relate $ and ^ is equivalent to 
setting g — (see eqn.(|3T|)). which is inconsistent with 
with allowing a second free function to modify Newton's 
constant. Making the choice g = uses up one degree of 
freedom, leaving us only a single free function with which 
to describe the system. 

The above reasoning is set within the confines of a 
second-order theory. We will now show that using a 
single function to relate $ and 'I' whilst maintaining 
Geff T^ Gq is equivalent to invoking a higher-derivative 
theory of gravity. To do this, let us consider the form 
that the tensor (Jt/^^ would take in a third-order theory. 
Its components would be: 



Ua 
Ue 
Up 

Ua 



Bok^ + Bi^ + lokt 

Cofc2$ + Cifc$ -1- 6*2$ + Jofc^f + Jifcf 
Di -Do- ~ Ki^ 



(35) 



The constraint equations for this system are given in 
Table [TTl We will continue to define the combination that 
modifies Go as g = —0.5{Ao + 3'HkBo), but note that this 
is no longer equal to HtDi as it was in the second-order 
case. Consider the case where we set Di = D2 = Kq = 
Ki — 0, that is, we use a single function to relate $ 
and ^. Through linear combinations of the constraints 





Origin 


Constraint equation 


1 


[Bl] $ 


io + HAo + kBo + nCo = 


2 


[Bl] $ 


ii + -Hyli + kAo + kBi + HCi = 


3 


[Bl] $ 


kAi + HC2 = 


4 


[Bl] f 


Fo + nFo + klo + nJo = 


5 


[Bl] f 


kFo + nJi=Q 


6 


[B2] $ 


Bo + 2HBo - ifcCo + |fc-Do = 


7 


[B2] $ 


Bl + kBo + 2-HBi - ifeCi + 2fcDi = 


= 


8 


[B2] -fe 


Bl - iC2 + 27^2 = 


9 


[B2] f 


Jo + 2-H/o - ifcJo + SfcA-Q = 


10 


[B2] f 


Jo-iJi + iifi=0 



TABLE II: Table of the constraint equations for the third- 
order metric theory specified in ijIVCI These can be generated 
using the formulae in Appendix [B] 



in Table In] we derive the expressions: 



Fo- 
5 = 



- SHkBi = 
iUklo = 



-Bl 



nj2 
tti. 



2k 



:Bi 



3E 



(36) 

(37) 

+ fc" ) (38) 



The first two of these expressions are the combinations 
that appear when we form the Poisson equation. They 
indicate that the potential additive modifications pro- 
portional to $ and r disappear; the format of eon. ([2^ is 
retained. Eqn. ([55|) shows that we can only have a modi- 
fication to the effective gravitational constant if i?i 7^ 0, 
and so from eqn. (|36p Ai ^ also. Using the third equa- 
tion in Table ITU C2 7^ in this case. Hence we are forced 

to include $ terms in U^ and Uq, and a $ term in Up. 
Since $ contains a first-order time derivative already (see 
eqn. p^ . the J/a in egn. pTj) will result in field equations 
containing third-order time derivatives - a higher-order 
gravitational theory. 

This result is a direct consequence of choosing Di — 
in C/s- Removing this constraint changes eqn. (j55)) to: 
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1 

2fc2 



Bl 



3E 



UkDi 



(39) 



which permits Bi — Q,g ^ Q, as we had in i jIVBI 

The above findings make sense within the context of 
the Lovelock-Grigore theorem [35l.l36|. which states that 
under the assumptions of four-dimensional Riemannian 
geometry and no additional fields, the Einstein-Hilbert 
action (plus a cosmological constant) is the only possible 
action that leads to local second-order field equations. 
In eqn. (|27l) the presence of Di/k in [/s means that this 
parameterization implies a non-local theory. This is not 
in itself problematic - nonlocal theories can arise when 
a degree of freedom has been integrated out, or elimi- 
nated from the action using an integral solution of the 
corresponding equation of motion. If Di = D^ = Ki = 



in eqns. psp there are no nonlocal terms present in the 
gravitational field equations, so we should not be sur- 
prised that the Lovelock-Grigore theorem prevents us 
from obtaining a second-order theory. In WII we will 
meet theories which evade the Lovelock-Grigore theorem 
in a number of different ways: by introducing new d.o.f. 
(e.g. scalar-tensor theory), higher-order field equations 
(/(i?) gravity), or through nonlocality and extra dimen- 
sions (DGP). 



D. Cases with 'XY' backgrounds 

The previous examples have all assumed that the 
background field equations are those of a Friedmann- 
Robertson- Walker metric plus standard cosmological flu- 
ids. We now relax this assumption and consider theories 
which modify the Einstein field equations at both the 
background and perturbative levels. It is well-known that 
any modification to background-level field equations is in- 
distinguishable from the effects of a dark fluid [37] ; hence 
we can write any background equations as the standard 
FRW ones with an additional energy density and pres- 
sure, see eqns.® and ^. We will refer to such theories 
as having 'XY backgrounds'. 

Any extension to GR must preserve the property of dif- 
feomorphism invariance. Invariance under passive diffeo- 
morphisms corresponds to the familiar principle of gen- 
eral covariance. Applying a passive diffeomorphism will 
generally result in field equations which look different 
to those in the old co-ordinate system. In contrast, in- 
variance under active diffeomorphisms requires that the 
actual form of field equations remains unchanged by a 
gauge transformation. In Table IIIII we list the gauge 
transformations for relevant variables. In practical terms, 
gauge form-invariance means that the extra terms that 
appear under a gauge transformation must cancel each 
other (using identities from the zero-order field equations 
if need be) . This places tight restrictions on our form for 

To see how this happens in ordinary GR, consider the 
linearly perturbed '00' component of the Einstein equa- 
tions (that is, eqn. p^ with C/a set to zero). When 
we apply a gauge transformation the left-hand side ac- 
quires a term ^SHES^/a. This is cancelled by the trans- 
formation of S on the right-hand side, provided that 
E = Ep -f Ep = dt-KGa? Ylii Pi(l + ^0' i-6- provided that 
the zeroth-order equations are satisfied. This is why 
we were only able to use gauge-invariant potentials in 
6U^u in CVTI and fJTVBJ if we don't aher the zeroth- 
order equations, adding anything else breaks gauge form- 
invariance. 

Now that we wish to consider XY backgrounds 
this procedure no longer works, because E ^ 
S-kGqo? Yli Pj(1 + Wj)- We must add a new term to C/a 
that will produce a part like -T-La^ {X -\-Y)£^ under a gauge 
transformation. Only then will the gauge-variant parts 
cancel by virtue of the zeroth-order equation. 
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Metric variables 


Fluid variables 


n. 


5^S-^{l + w)H^ 


^^^ + ^\i + H^-M 


e^e + ^^ 


^^/3+i[6-«C-2fcVl 


n -^ n + ^[w - SwHii + w)]^ 


u -^ u + ^i) 


s^ s 


v^v + U 






Components of SG^iu 


EA^EA-'^H{EF + Eii)£, 


Ee^Ee + ^ {Ep + Er) ^ 


Ep^Ep+'^{ER-2nER) 


Es -^ Et. 



TABLE III: Behaviour of metric and fluid variables under 
infinitesimal diffeomorphisms generated by the vector field 
Cm = fl(— C)^i^)- Note that the shear E is gauge- invariant. 



As a toy example, consider a simple theory which 
modifies the zeroth-order equations solely by introduc- 
ing time-dependence to Newton's gravitational constant. 
Following the notation of previous sections, we can write 
the sum of the modified Fricdmann and Raychaudhuri 
equations as: 



E 



S7^^-^^^^a^Yp^il + w^) 



1 - gb{a) 



/ , I 



(40) 



Rewriting this in the form of ordinary GR (and hereafter 
suppressing the argument of gt): 



E = SnGoa^ ^ /3<(1 + w^) + gbE 



(41) 



explicitly first-order in derivatives, resulting in second- 
order equations. 

Uq is treated analogously to Ua, and has the form: 



Ue = 2gb 



kr + \ev 

4 



Bnfc$ 



(43) 



Combining egns. p^ and (H5|) . we find that the Poisson 
equation has the same form as it did in the case with un- 
modified background equations, eqn.([29]). We define the 
combination g = —1/2 (Ao -t- 3HkBo) as we did in i jIVBI 
We will assume that the modifications to Newton's grav- 
itational constant appearing in the zeroth-order and per- 
turbed equations are the same, i.e. gi, = g, noting that 
we have not formally proved this to be the case. 

Once we have deduced the form of Ua and Uq , Up and 
Us can be found using the Bianchi identities (eqns. fTTj) 
and (|12[) ^ . Eliminating the free function ^o in favour of 
g and Bq, these are: 



HUp = 



3$ 


kUB, 


+ kB^L + 


5^MS 


+3^ 


(E k^\ 


+6 {tn + gL) 


kV^-EV 
4 


^'-K. 


kV^-EV 
4 





3 






-E~gnV 



H~g) 



(44) 



from which we can identify a'^{X+Y) — g^E (see eqns.© 
and ([5])). A possible form for Ua is then: 



Ua = -QgbT-L 



kV 



-EV 



Aok^^ 



(42) 



Unlike the second-order example of ^IVB[ Ua now con- 
tains a r term. The offending second-order derivative 
of the scale factor present in T is cancelled by the term 
proportional to EV. We did not have the freedom to 
add such a term in the case of unmodified background 
equations, because V is gauge- variant. 

Using ean. (l42t in egn. dTOl) and the transformations 
given in Table Hill it can be verified that the gauge- 
variant parts cancel by satisfying ean.(|i(71). Note that 
the need to have gauge-invariant second-order equations 
has totally fixed the first term in egn. lH^ : all freedom re- 
sides in the gauge-invariant part of Ua via the function 
Ao. 

As a quick sanity check, one can verify that the field 
equations remain second order in the conformal Newto- 
nian gauge. In this gauge V = V = 0, E ~ —"^ and 
/3 — —6$, so $ and 'P reduce to their familiar coun- 
terparts $ and VP that appear in the linearly perturbed 

FRW metric (recah that f = l/fc($ -I- -H*)). Ua is then 



mk 



Us 



-— ffcr-$ 

where L^U^ +% 



2k{g + n~g)-iB„[- + - 



E 



fc2 

3 ' 



(45) 



As expected, we find that Uy, contains only gauge- 
invariant perturbation variables. This must be the case 
since all other terms in ean. (j22p are gauge-invariant, so 
there is nothing to cancel against. 

If we express the last term in Us in terms of ^ we 
find that we can write the relationship between the two 
potentials as $ = f?siip(a:fc)^ for this toy example. We 
have already chosen our two PPF functions to be g and 
Bq, so 77siip is simply a particular combination of these: 



??slip 



H(l-g) 



H(l-.g)-,^ + lfi(f + f) 



(46) 



The key result of this toy example is that (for purely 
metric theories) SUf^i, can have a more complex form 
when the background equations are not standard FRW 
plus standard cosmological components (baryons, CDM, 
etc.) The hierarchy of constraint equations becomes more 



complex due to the non-zero X and Y terms in eans. pTj ) 
and p^ . We can use the principles of energy conserva- 
tion, gauge-invariance and second-order field equations as 
a shortcut to the correct forms; the same results would be 
obtained by solving the hierarchy of constraint equations 
directly. 



V. THE GENERAL PARAMETERIZATION- 
EXTRA FIELDS 

The formalism we have developed so far is only ap- 
plicable to purely metric theories. Yet, as discussed in 
jjlll the majority of modified gravity theories introduce 
new degrees of freedom, often as additional scalar, vec- 
tor or tensor fields. It is not immediately obvious that 
the behaviour of these theories can be be encapsulated 
by a either a (Q,?7siip) or (5,C) parameterization. There 
is a risk that we might develop a model-independent for- 
malism that does not map onto most of our well-studied 
theories. 

In WII we will study several example cases, chosen to 
be representative of common classes of modified gravity 
models, and ask whether they can be expressed in the 
two-function format of quasistatic PPF. Attempting to 
map disparate theories onto a single framework is only 
plausible if those theories share some common features. 
Hence, before turning to specific examples, we wish to 
consider what general statements can be made about the 
structure of the the field equations in theories with extra 
degrees of freedom. 

Let the scalar perturbations to an extra degree of free- 
dom be denoted by x-, e.g. if the new degree of freedom is 
a scalar field </> then x = 5(j). (We have succumbed to the 
common but unwise choice of terminology by referring 
to 'scalar perturbations', even though the new degree of 
freedom may itself be a vector or tensor field. 'Spin-0 per- 
turbations' would be a better choice of terminology [38|). 
If we are to obtain second-order field equations then we 
know that C/a and Uq can only contain the perturbations 
X and X- 

More generally we can introduce multiple new d.o.f. 
and denote their scalar perturbations by the vector Xi 
with components x ■ The perturbed field equations in 
a general gauge are awkward and rarely used; hence we 
will specialise to the conformal Newtonian gauge for the 
remainder of this paper. The relevant expressions for 
scalar-tensor theory ( WI A[) are presented in a general 
gauge in Appendix VK\ 

In the conformal Newtonian gauge the 00-element of 
the tensor 5U fj_i, can be represented as: 

[/a = k'^a^x + kdjjl+ Aok^<^ + Aiki^ 

+FokH + Fife* + ^5p + ^Sp (47) 

where 5p is the total energy density fluctuation of stan- 
dard cosmological fluids (similarly for 0, 11 and S to be 
used shortly). q?o and cJi denote vectors of functions 



with the same dimensionality as x- <5i, Aj, Fi and Mi axe 
functions of background quantities such as p and a; these 
dependencies have been suppressed for clarity. Note that 
ean. (|T7)) has the form indicated schematically in eqn.®. 
The term MiSp represents a modification which depends 
on the rate of change of the density fiuctuations of or- 
dinary matter. Whilst formally this term is permitted 
to be present in C/a, we are unaware of any theory of 
modified gravity that results in a perturbed 00-equation 
with a term like this. Theories employing a chameleon 
mechanism introduce modifications to GR that depend 
upon the environmental matter density, but not upon its 
rate of change . We will therefore choose Mi = in what 
follows. The 5p term in eqn. (|47|) can be eliminated using 
the relation: 



SirGoa'^dp — E/\ — Ua 



(48) 



Ea can be expressed in terms of $, $ and ^ in the con- 
formal Newtonian gauge, see egns. fTS) ). Rearranging and 
redefining the coefficient functions, we obtain: 



C/a 



k'^^ox + kdjj( + Aok'^<i> 



Aifc$ 
+ Fok'^-9 + Fikii (49) 

This procedure has enabled us to eliminate energy 
density fluctuations of ordinary matter from C/a. If the 
new d.o.f. are not coupled to ordinary matter then Sp 
does not appear in C/a anyway. It is worth reminding 
the reader that in this section we are working in the 
conformal Newtonian gauge, so <I> and $ should not 
be confused with their gauge-invariant counterparts <& 
and ^, which already contain first- and second-order 
derivatives respectively. If we were to consider eqn. (j49p 

in a general gauge we would find that <& and ^ only 
appear in the combination F. The ^ in ean. (|49p would 
become — S (recall that in the conformal Newtonian 
gauge $ = $ = — ^x and ^P = ^I* = — S, using egns. p^ 
and pi|) ). 

We apply a similar treatment to the remaining compo- 
nents of SUni, . Recall that Up and C/s are permitted to 
contain second-order terms such as $ and x- Terms such 
as 0,n,n, S and S are discarded to maintain contain 
consistency with our treatment of C/a; we stress again 
that this is done only on an intuitive basis. The result- 
ing expressions for SUf^i, (together with eqn. pg)) ) are as 
follows: 

Ue = /cAfx + /3?"f +So/c$ + Bi$ + /ofc* + /i^ (50) 



Up = k'%^x + k^iX + l2X + Cok^^ + Cik^ 



-C2$ + Jo^ * + -^1^* + ■^2* 



^, 



1 



jr,^ 



1 



Di 






(51) 



(52) 



where /3i, . . . Si and Bi, . . . Ki denote functions of back- 
ground quantities. 
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The Bianchi identities then give us two constraint 
equations couphng terms in x i^i^ ^^nd their deriva- 
tives. In contrast to the previous section we can no longer 
set the coefficients of the each term to zero individually. 
In the case without extra fields this was possible because 
all our variables were non-dynamical, so obtaining evo- 
lution equations for them would be unphysical. But now 
that extra fields appear in SU^,^, the Bianchi identities 
yield equations describing how the metric variables re- 
spond to the set of perturbations x- Therefore we no 
longer have a hierarchy of constraint equations for the 
cocffcicnts ao, . . . K2 that allow us to reduce them down 
to two functions. This is not problematic in itself. To 
map a specific theory onto the parameterization we can 
simply pull the necessary coefficients out of the perturbed 
field equations. We will see shortly ( WI|) that in many 
cases that these are relatively simple functions. 

To constrain a general parameterized theory using 
Markov Chain Monte Carlo (MCMC) analysis we instead 
choose some sensible ansatz for the functions ao, . . . -ft'2. 
For exarnple. a Taylor series up to cubic order in f^A was 
used in [3^, |40| ; the MCMC then constrains the coeffi- 
cients of the terms in the Taylor series. Rigidly fixing 
the format of the parameterization in this way means 
that we simply have to constrain real numbers. This 
simplicity is a key advantage of explicitly parameterizing 
for the new fields as in egns. ipMS^ . The alternative ap- 
proach - absorbing the new fields into an evolving Gcs 
and slip parameter where possible - will give GcS and C 
very complicated forms that are difficult to parameter- 
ize (for example, see eons. ([751) and (1^51) ). The trade-off 
is that our method requires considerably more than two 
coefficient functions. We expect that some of these will 
be well-constrained by the data, others less so. 

In the case of just one or two new d.o.f., the system 
consisting of the Einstein equations, the two conservation 
equations for ordinary matter and two Bianchi identities 
for the [7-tensor can be solved. In order to avoid a con- 
tradiction, the Bianchi identities for the [/-tensor must 
be equivalent to the equations of motion for the extra 
degrees of freedom (obtained by varying the action with 
respect to the extra fields or similar). Futhermore, when 
only a single d.o.f. is present the solutions of the two 
components of the Bianchi identity must be consistent 
with each other. 

When more than two new d.o.f. are present the Bianchi 
identities do not provide sufficient information to solve 
the system, and we must supply additional relations be- 
tween the new d.o.f, metric variables and matter vari- 
ables. With our goal of an abstract, unified framework 
in mind, we will introduce a general structure to tackle 
such cases. We make the conjecture that one can use the 
'generalized dark matter' (GDM) formalism developed by 
Hu |41| in order to obtain the necessary closure relations. 
GDM is a phenomenological model in which specification 
of three parameters - an equation of state, a rest-frame 
sound speed and a viscous sound speed - suffice to recon- 
struct the full perturbed stress-energy tensor of a fiuid. 



Cold dark matter, radiation, massive neutrinos, WIMPs, 
scalar fields and a cosmological constant can all be re- 
covered as limiting cases of GDM. In our case the d.o.f. 
parameterized as GDM may be genuine fiuid components 
(e.g. scalar or vector fields), or effective fiuids (eg. the 
scalaron of f{R) gravity, the Weyl fiuid of DGP gravity 
- see WI|) . For example, in eqn.(j3ni) we identify the extra 
d.o.f. with an energy density perturbation: 



d^X + "I'x = SnGoa'^PE^E 



(53) 



where the 'hat' symbol indicates that we have folded 
in the necessary metric perturbations to make gauge- 
invariant versions of x ^-iid %. Similarly we can con- 
struct gauge-invariant x, x ^nd x from the terms in 
cans. ([50)1 - ([5^ . and identify these with velocity, pressure 
and anisotropic stress perturbations respectively. A sub- 
script E will be used to indicate these effective perturba- 
tions. 

The GDM formalism then provides a way of reducing 
these four fiuid perturbations to just two, which can then 
be related to the metric potentials via the perturbed con- 
servation equations. These are |41| : 

weSe 



1$ 



..]$ 



-(1 + w'b)( 

-[...]4>-i-[. 



k^0f 



3$ 



1 



1$ 



* 



We 
. .1* 



SHweTe + 
[...]* (54) 



-n 1 



3c' 



ad; 



cIJe + weTe 2, 



1$ 



..]$ 



(1- 
.1* 



we) 



'^ + [...]^ (55) 



where the square brackets denote combinations of the 
functions Bi, . . .Ki. We should remember that these are 

really second-order equations, due to the x in ean. (l53|) . 
In can. ([55]) the pressure perturbation of the effective 
fiuid has been decomposed into an adiabatic and a non- 
adiabatic part: 



IIb = c^d^B + weTe 



(56) 



We have adopted common notation by using V e to rep- 
resent the dimensionless non-adiabatic pressure pertur- 
bation; this should not be confused with our metric po- 
tential F. The adiabatic sound speed is fully determined 
by the equation of state parameter we'- 



1 



-"ad 



We 



We 



in I + WE 



(57) 



The non-adiabatic pressure perturbation is specified by 
introducing a parameter c^^, interpreted as the sound 
speed of the fiuid in its rest frame: 



we^e = 



■-off 



cl^){5E + m{i + WE)eE 



(58) 



A third and final parameter is needed to relate the 
anisotropic stress, S^;, to the velocity perturbations. 
This is the viscosity parameter c^ 



r.2 

"■vis" 



{l + WE)[t 



)(= 



3-HSf 



WE , 
WE 



4.cl:,k'0E 



(59) 
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By combining eans. ((57|) - (l59l) with cans. ([54) 1 and (|55|) we 
can obtain two equations relating any two of the fluid per- 
turbations to the two metric potentials (although we note 
that the presence of the E^; term in eqn. (j59p might make 
this step non-trivial). These expressions will contain the 
three GDM parameters {wE,clff,cl^^}. An equivalent 
three-parameter framework was studied in the context 
of dark energy in [42| , and constraints from current and 
future data sets were investigated in [4j] . 

We have already mentioned the degeneracy bewteen 
modifications to gravity and fluid components in the 
zeroth-order field equations [33] ■ An explicit example 
of this is presented in [4J] for the case of Eddington- 
Born-Infeld gravity. We stress that the effective fiuids 
representing the extra d.o.f. at the background and per- 
turbed levels need not have the same properties. Indeed, 
it is most likely that they will be different. If they are 
not, we have no way of deciding whether dark energy is 
really a modification to gravity or a dark fiuid, even using 
observations that reflect the rate of growth of structure 

As described above, for theories with two or less extra 
d.o.f. the GDM prescription is not strictly needed; the 
system of equations is already closed. But if we wish to 
constrain modified gravity in a model-independent way 
we cannot make assumptions about the number of new 
d.o.f. introduced. The GDM approach allows us to ob- 
tain model-independent closure relations at the expense 
of introducing three new parameters, which would need 
to be constrained in a MCMC analysis. 



VI. EXAMPLES 

Thus far our discussion of the PPF parameterization 
has been purely formal. Studying some specific cases 
of gravitational theories should help to consolidate the 
ideas outlined in this paper. The four examples pre- 
sented here are chosen to represent some common classes 
of theories of modified gravity; a comprehensive review of 
many other theories is presented in [Ij. A simple scalar- 
tensor theory and Einstein- Aether theory represent the- 
ories that introduce additional fields to GR (scalar and 
vector fields respectively). Theories with additional ten- 
sor fields such as Eddington-Born-Infeld gravity [4Jj Eg] 
and bimetric theories [47[ also belong to this broad cate- 
gory. f{R) gravity is studied as an example of a higher- 
derivative theory; Hofava-Liftschitz gravity pSi - ISOf and 
Galileon theories jSlI, [52] also fall into this class. Our fi- 
nal example is Dvali-Gabadazze-Porrati gravity (DGP) 
|53l [SJ], which we study as a representative higher- 
dimensional theory. Whilst DGP itself is now disfavoured 
by observations, it incorporates features common to other 
braneworld theories [5a, [5g| such as Randall-Sundrum 
models [53, [sl] and cascading gravity [S^ . 

For each theory we will consider the extent to which 
the gravitational field equations can be modelled by the 
PPF parameterization. Is it possible to describe such a 



rich variety of theories using only two functions? The 
usual arena for PPF is the quasistatic limit, in which 
time derivatives of perturbations can be neglected rela- 
tive to spatial derivatives. This is the dominant regime 
for measures of the rate of structure growth, such as weak 
lensing and peculiar velocity surveys. However, we are 
also interested in using constraints from the Integrated 
Sachs- Wolfe effect, which requires consideration of scales 
above the quasistatic regime. There has been much work 
recently highlighting the importance of a correct rela- 
tivistic treatment of large scales, for both the theoretical 
and observed matter power spectrum |60l463l |. We will 
consider whether the PPF parameterization can be ex- 
tended to this regime. 



A. Scalar-tensor theory 

Scalar-tensor theories, first considered by Jordan, 
Brans and Dicke [sj , modifiy GR by introducing a scalar 
field which couples to the Ricci scalar in the gravita- 
tional action. The concept is closely linked to that of 
quintessence, in which a scalar field is used a dark en- 
ergy fiuid-type component but without the explicit non- 
minimal coupling to the Ricci scalar. Particle physics 
is in no short supply of candidate scalar fields, and the 
reduction of higher-dimensional theories to effective four- 
dimensional field theories also gives rise to candidate 
scalars (moduli). However, finding a field with exactly 
the right properties to account for dark energy has proved 
difficult. There is an obvious aesthetic appeal in connect- 
ing the scalar field with the infiaton; however, they need 
not a priori be the same field, and introducing such a 
connection creates further obstacles to constructing vi- 
able models. 

In general, three functions are required to specify a 
scalar-tensor theory: the coupling to the Ricci scalar, 
F(0), a potential for the scalar field, U {(/)), and a func- 
tion multiplying the kinetic term of the scalar, Z{(j)). 
However, it is possible to reduce F{(j)) and Z{(f>) to 
just one function through a field redefinition, result- 
ing in F{(l)) = <j) and Z{(f>) ^oj{(l))/(j) ^. The choice 
a;((/)) — constant is termed a Brans-Dicke theory, which 
recovers GR plus a cosmological constant in the limit 
ijj — >■ oo, U{(f)) — > A. Measurements from the Cassini 
spacecraft constrain cj > 40, 000 (2cr limits) in the solar 
system ^^. 

A generic property of scalar-tensor theories is that they 
result in a time-dependent gravitational 'constant'. This 
is precisely one of the features of the PPF formalism 
(through g in our parameterization), which gives hope 
that scalar-tensor gravity might be fully accommodated 
by PPF. For simplicity we will focus on a scalar-tensor 
theory with variable coupling w(0) and no potential, 
working in the conformal Newtonian gauge. The action 
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in the Jordan frame is: 



^=16^/^'^^ 



'^i?_^(V^)2' 



Sm[^P'^'^,g^.. 



(60) 



where Sra[ip''^\gpLij\ is the matter action and V'^*-' the mat- 
ter fields. Varying this action with respect to the metric 
yields the gravitational field equations |67| : 



^ ^u — 



SvrGo^ , Lo{cp) 



'^^- + ^ (v^0V.0 - ig^.(V0)^ 



-- (V^V^(/) - g^^ D ( 



(61) 



where D = g'^'^V^Vi/. Gfj,i, is the usual Einstein tensor 
of GR and the scalar field has been rescaled so that it 
is dimensionless, cj) — > (l)/Go. In a smooth, unperturbed 
FRW universe this gives us the background equations 
(c.f. equs.® and ©): 



Ef 



SttGq 



a^Y.p^ + U^)^-^'^: 



Er = — —a^ 2^ P, + -a; ((/))■ 



n- 



(62) 
(63) 



Variation of the action with respect to (j) gives the equa- 
tion of motion for the scalar field [63 : 



U<t) 



2w((/)) -I- 3 



87rGoaX^-^^(V0)2) (64) 



Rewriting eqn. (pT|) in the form of eqn.Q indicates that 
the form of C/^jy must be: 



U^, = G^,(l - 0) + ^ (v^0V,0 - i<?^.(V0)2 



-V^Vy0- gpi/D' 



(65) 



Linearly perturbing this expression (and raising an in- 



dex) will give us U^ ,Uq,Up and C/s . We obtain: 



Ua = Ea{1 -(/)) + ^[3(1)] + ^ 

1 duj{(j)) (/)2 1 



&n4> - ijj{(i>)^ 



^(0)- 



3-H' 



-in + uj{(j))^ 



Ue = ^e(l 



a;(0)-j--H 



Up = £;p(i -(/)) + $ -60 



* 



-60-6H0-3w(0)-- 



-L.(0) 



02 3duj{(t>) 02 



(66) 



50 + 1'[-0] (67) 



^ [-30] 



3-^2 + 6-H + 2fc2 



3H + 3w(0)^ 



C/s = ^s(l- 



(68) 
(69) 



Specialising to the conformal Newtonian gauge (in which 
V=0) means that we lose the time derivatives in eqns.(fT5|) 
and ([TB)) . Hence the appearance of ^ in Ua and Uq 
above gives us no cause for concern; the Bianchi identi- 
ties will remain second-order equations. If we had kept to 
a general gauge additional terms in V, e and v would be 
present in the above expressions, ensuring that ^ only ap- 
peared within the combination $ + H\t' and that any a/ a 
terms cancel. The (lengthy) corresponding expressions 
for a general gauge are displayed in Appendix [Xj There 
we also demonstrate that the general-gauge expressions 
still obey the constraints of yielding second-order, gauge- 
invariant perturbed field equations. 

We can see that ean. (l55)) has the form indicated in 
ean. (l49l) : in this particular case we can pick out the co- 
efficient functions: 



Ao = -2(1- 
30 



Ai 



6Hfe(l 



Mk 



Fo = -6^2(1-0) +6^0 
Fi = 



C^(0) 02 



ao = 



k '' k?' (j) 
1 dw(0) 02 1 -^2 



2fc2 



^-7Tr2^('^)72-l-3^A 



2fc2'"^'^^02 ^ ""-fc 



ai = —SHk + 



^(0) 



(70) 



Similarly one can read off the coefficients Pi, . . .£i and 
Bi, . . . Ki for scalar-tensor theory by matching eans. lfSO) )- 
(15^ with eans. ((S7)) - (|M|) . 
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Using the expressions for Ua and Uq in eqns.® we 
can form the modified (Fourier-space) Poisson equation. 
Some of the terms in 6(j) and (50 can be combined to form 
the gauge- invariant density perturbation of a scalar field: 

p^A^ = 6p^ + in (p0 + p^) 00 = ((50 -4)^ + in (5(/)) 

(71) 



The Poisson equation for scalar-tensor gravity is then: 



+3- 



$ + -H* 



1 (/)2 /w((?i) dw((/)) 



fc^+6H^ + - — 



Ai is a gauge-invariant density perturbation to matter 
(CDM, baryons, radiation), and we have redefined the 
scalar field so as to pull a factor of i-KGoo? out of the 
second term. We have chosen to write the Poisson equa- 
tion in this form because it delineates the extra terms 
that arise in a scalar-tensor theory compared to uncou- 
pled quintessence. We reach a quintessence-like limit by 
setting ijj{<f) = (j) and removing the 4>-R coupling in the 
action. As a result, the last two lines of ean. ([7^ do not 
arise and the prefactor is just SttGoo^ (because we don't 
rescale cj) — >■ (j)/Go). It is the last two lines that offer the 
distinction between quintessence and scalar-tensor grav- 
ity. The slip relation is: 



S6 
$-* = — 



(73) 



It is useful to consider some simplifying limits of the 
above expressions. The 'smooth' limit, in which the per- 
turbations of the scalar field are negligible relative to the 
matter perturbations, gives us: 



-2e<p = 8TT^a^yp,A, 



$ + -H* 



(74) 



The second term on the right-hand side cannot be ab- 
sorbed into the first term without giving GcS an unde- 
sired environmental dependence. 

However, if we take the quasistatic limit of eqn.(|72l) - 
that is, we neglect time derivatives of perturbation vari- 
ables and take k ^ n -we find that the modifications to 
the Poisson equation can indeed be repackaged as a mod- 
ified gravitational constant. In this limit it is possible to 
write the relation between the two metric potentials as 
^ = 'ysiip^; e.g. for the choice w((^) — cj) the form of ?7siip 
is particularly simple [6a, l63 : 



1 



?7siip 



(75) 



Combining cans. ([72) 1 and (l73l) with the slip relation, we 
obtain: 



^k^^^S^^a^lr^p^A, 



1 / 1 

?/siipy 2fc2^ 



uj{cp) 



bA^ X (76) 



Vslip 



uj{(j)) duj{(j)) 



■ -i -1 



in 

fc^?ysiip 4> 



where the expression in square brackets gives the time- 
and scale-dependence of Newton's gravitational constant. 
As mentioned earlier, it is important that we also con- 
sider the (super)horizon-scale limit of these theories for 
correct treatment of their predicted effects on the matter 
power spectrum and large-angle CMB power spectrum. 
Taking the limit k ^n allows us to neglect the k'^S<j>/(f> 
term of eqn. (j72p . but there are no other obvious simpli- 
fications unless (5(/) is negligible on these scales. It seems 
that in order to cope with the superhorizon limit we need 
a parameterization that allows for additive modifications 
to the Poisson equation, as well as a modified GcS- Of 
course, there is no barrier to using the standard PPF 
format to model part of the modifications, but we need 
to remember that the correspondence between parame- 
terization and theory would no longer be exact on large 
scales. 



B. f{R) gravity 

Another commonly-studied theory is f{R) gravity, for 
which the action is a general function of the Ricci scalar, 
/(i?); see [TO] for a detailed review. Two approaches 
to f{R) gravity are possible. In the metric formulation 
the affine connection T^ is defined in terms of the met- 
ric components in the usual way, and gravitational field 
equations are obtained by varying the metric with respect 
to g^^ only. In the Palatini formulation of f{R) gravity 
the connection and the metric are treated as independent 
variables and the action is varied with respect them indi- 
vidually. One finds that in the metric formulation there is 
a propagating scalar degree of freedom /^ = df{R)/dR. 
Since the Ricci scalar is constructed from second deriva- 
tives of the metric, the kinetic term of the scalar degree 
of freedom contains fourth-order derivatives and hence 
metric f{R) gravity corresponds to a higher-order the- 
ory. 

The presence of a scalar degree of freedom within f{R) 
theories can be made explicit by applying a conformal 
transformation that maps f{R) gravity onto a scalar- 
tensor theory. Metric and Palatini f{R) gravity map 
onto Brans-Dicke theories with oj = 0, — | respectively 
[70|. The scalar field arising under the conformal trans- 
formation is sometimes referred to as the 'scalaron'. In 
the conformally-transformed frame (the Einstein frame) 
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the scalaron acquires a coupling to the matter fields V'*-*\ 
leading to non-standard conservation equations. Hence 
the non-transformed frame (the Jordan frame) is re- 
garded as the physical frame in which observations are 
made. 

The action in the Jordan frame is; 



S = 



1 



WttGo 



d'xy^/(i?) + 5,„[V^^\.9p.] 



(77) 



Variation with respect to the metric leads to the field 
equations: 

(78) 

This can be rewritten in the form of eqn.(|4|) with the 
following expression for f/^i^: 



above expression 



Ua = ^A(i-/fl,) + $[3/fl] + * 6ni 





+X 


sn-P 


+ 


X [-m 






(80) 


Ue 


= Eeil-fR)+x[-H]+x + ^[-fR] (81) 


Up 


= Ep{l-fp) + <i> 


-6fR 


+ * 


-6/fl - 6nfR 


(82) 




+^ h3/fl] + X 


-3H - 6H^ + 2fc2 


+ x[3H]+3x 


C/s 


= E^{ 


!-//?) + 


X 










(83) 



The similarity between the sets of eqns. (l66l) - (|69)) and 
(|80| -([83 |) is immediately apparent, suggesting the identi- 
fication of fp as the scalaron. However, in f{R) gravity 
the 'extra' d.o.f, x, can be expressed in terms of metric 
potentials, x ^nd x ^re given by the expressions: 



U^. = i?,,.(l-/fl)-25M-(-R-/(-R))+V„V„/i,-.g^,n/fl 

(79) 

We will let x denote the perturbation to the extra de- 
gree of freedom, ie. x = ^{Ir) = Irr^R- Perturbing the 



J 



f3RRSR + fRR{SR) 



(84) 



X 



6R (UrR^ + hRR) + 2hRR{5R) + fRR{6R) (85) 



where fm indicates the fourth derivative of / with re- 
spect to R etc., and 



a'^SR 



a'iSR) 
a\SR) 



-4fc2$ - 18H$ - 6$ - 2 {6n + en^ -k^j^^- 6-H^ (86) 

-2a'^n5R - $ Uk^ + isn) - mn^ - 6$^^) - 12* (n + 2%%^ - * (im. + un^ - 2kA - en^ (87) 

-SRa^ {m + AuA - Aa^l-i{5R) - 18H$ - $ {iQlk. + AkA - l^U^^^^ - 6$^'') - 12* {u^^'i + 2%^ + 2%%] 
-* (30?^ + AmH\ - * hm + 12-H^ - 2kA - 6-H*(3' (88) 

I 



The coefficients corresponding to the general expressions 
in ^V] can be computed. For example, matching onto 
eqn.dUl): 



Ao 



,/. 



RR. 



2 (1 - Ir) - 4^ [3n-k' + 6n 



Fo 



-ml{i-fR) + Q^JR + &^-^ 



un^nl 



-enln - ml + enkHk + m^ + m- — 



+6^nR 



6nt + 6 



-Hk 



-1 







F, = 18^ 


1 


= 6Hfe(l /i,)+3^/'+6^^^J^'=(5fc2 ISn') 
K, a 


F2 = 18H' 




+54p^n^R 

ka^ 


4^^ = «« 


2 


= 6^ (k^ + sn^ - 3h) + is^nii 




3 




The slip rela1 



18^ (2/flfiH + hR-HR 



2JRR 



a" 



„(*) 



(89) 



The slip relation and modified Poisson equation for f{R) 
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gravity become: 

^ - -^ = 2L 
Ir 

JR ; JR 



(90) 



$ + H* 






-E + k- 



(91) 



Again we see that these expressions are largely similar 
to the scalar-tensor case with a;(0) = 0, fn replacing 
(j) and X = ^4'- However, we note that the similarities 
are aesthetic only; the equivalence between scalar-tensor 
theory and /(i?) gravity is only formally realised under 
a conformal transformation as described above. Scf) rep- 
resents a perturbation to a new field that is genuinely 
additional to GR, whereas x is really a placeholder for 
the combination of metric perturbations in eqn. (|86p . 

In the 'smooth' limit % — >■ we obtain the analogy of 
eqn.dlll): 



_2fc2$^8^5^a2^ 



Jr 



Pi^i 



Jr 
Tr 



$-fH* 



(92) 



with the same second background-dependent term as the 
scalar-tensor case. 

For measures of late-time structure growth we are pre- 
dominantly interested in the quasistatic regime. Using 
X = Irr^R a-nd eqn. ()86p we can replace x in the Pois- 
son equation and slip relation by sequence of the metric 
potentials and their derivatives, just as we laid out in 
eans. (P5) )- (P^ . For example, in the quasistatic limit the 
slip relation becomes: 



$-* 



fRR 

fRa' 



(.(- 



6H -I- fc2)* - 4fc2$ 



(93) 



This matches the result of Pogosian & Silvcstri 71] if 
we further neglect the T-l term. They expressed the slip 
parameter of f{R) gravity in the quasistatic limit as: 



?7slip = 



3 + 2Q 
3 + 4Q 



(94) 



where Q is defined as the squared ratio of the Compton 
wavelength of the scalaron to the physical wavelength of 
a mode: 



Q 



, fc^ Jrr 
'a^ Jr 



Ac 
A 



(95) 



In this limit the Poisson equation can be rewritten as: 



Go 



fc2$ ^ STT^a^^p^A 



Ir 



(96) 




IrU 
iRk^ 



3 E 

2fc2 



1 

'7slip 



The term in square brackets controls the time- and scale- 
dependence of the effective Newton's constant. 

It is no surprise to find that f{R) gravity behaves sim- 
ilar to scalar-tensor theory in the superhorizon limit. We 
can neglect k'^ in ean. ([9T|) and eliminate x using ean.(|90|. 
but the resulting expression will still have terms in ^ and 
$ on the right-hand side that can not be written as a 
Gcff. The conclusion is analogous to that of iJVI Al the 
mapping between a (Geff,7ysiip) parameterization is ex- 
act in the quasistatic limit, but only approximate in the 
large-scale regime. 



C. Einstein- Aether theory 

The emergence of string theory as a candidate theory 
of quantum gravity leads to the possibility that space- 
time coordinates are non-commutative [72] . Under these 
circumstances Lorentz symmetry may be violated. Hence 
much effort has been invested in exploring ways in which 
Lorentz violation can be implemented without marring 
the key successful features of GR, such as general covari- 
ance. A minimal way to do this is to introduce a vector 
field into the action, which defines a preferred reference 
frame at every point in spacetime. In Minkowski space- 
time one can simply introduce a constant vector field 
into the background, but in a curved spacetime this is 
not possible; the Lorentz-violating vector field must be 
promoted to a dynamical field derived from a generally 
covariant action. If invariance under three-dimensional 
spatial rotations is preserved then the vector field must 
be time-like. Einstein- Aether theories [73| introduce such 
a vector field ('the aether') of unit length, which is con- 
strained not to vanish so that Lorentz violation is main- 
tained even in a vacuum. The unit, time-like nature of 
the vector field is enforced through means of a Lagrange 
multiplier in the action. 

TeVeS |74| is a well-known example of another theory 
that contains Lorentz-violating vector fields. When for- 
mulated in the Einstein frame, TeVeS introduces both a 
unit, timelike vector field and a scalar field, and employs 
a free function to ensure that it reduces to Modified New- 
tonian Dynamics (MOND) [Tal in the non-relativisitc 
limit. Upon transformation to the Jordan frame - in 
which particles follow geodesies of the metric - the unit 
length of the vector field is not preserved. It was shown 
by Zlosnik [73] that the scalar field present in the Ein- 
stein frame is absorbed by the vector field in the Jordan 
frame, and dynamically determinines the modulus of the 
vector. Hence TeVeS is equivalent to an Einstein- Aether 
theory in which the aether field has variable length. Our 
study of the compatibility of Einstein- Aether theory with 
the PPF framework therefore has implications for TeVeS 
also. 

The most general action for Einstein- Aether theories 
is: 



5 = 



^ j d^xT^ [R + Cea{9^\ A^)] 



+ Sraiw''' ,g^lu 
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(97) field equations are: 



where the aether Lagrangian is 

(^EA {g'^^A^) = M^TilC) + \{A^A^ + 1) (98) 

M has the dimensions of mass. /C is a scalar formed from 
the kinetic term of the vector field: 



/C = Af~2^^^y^^7V^A'^ 



'^70- ~ 



Cl9"'^9jo 



C2<5"^^+C3^X 



(99) 



where Ci are dimensionless constants. Barbero & 
Villasefior [73] have identified special choices of a for 
which Einstein- Aether theory becomes equivalent to GR 
under a field redefinition. More generally, the Ci must 
obey certain restrictions if the linearized field equations 
are to be hyperbolic - hence admitting a well-posed ini- 
tial value problem - and exclude superluminal propaga- 
tion of aether perturbations and gravitational waves [73| ■ 
In general there is potentially a fourth term in ean.(|99p. 
but in the case of purely spin-0 perturbations this can be 
absorbed by suitable redefinitions of ci and C3. We will 
assume there is no direct coupling between matter and 
the aether; they interact only gravitationally. 

The gravitational field equations obtained by varying 
the Einstein- Aether action with respect to the metric can 
be written in the form of eqn.Q, where the modification 
tensor [/^^ is effectively the stress-energy tensor of the 
aether: 



U^. - 



V. [Tk (J(/^.)) - J\^A,) - J(^u)A'' 



'^kY^^u + -izQtiu^ + >^A^A^ 



(100) 



where round brackets around subscripts denotes sym- 
metrization with weight 1/2. We use the notation J-k = 
dT jdK,, and the following definitions: 

J\ - ^t^0A^ (101) 

V - ci (V"A^V„^, - V^A"V,A„) (102) 

Varying the action with respect to the aether field gives 
the equation of motion: 



V^ (J-K J^ J = 1\A, 



(103) 



Finally, varying the action with respect to the Lagrange 
multiplier A gives the constraint A^ A^ = — 1. The re- 
quirement of a spatially isotropic background fixes A^ = 
(1, 0, 0, 0). Using this in ean. (|103p enables us to solve for 
the Lagrange multiplier: 



A = --A^V. {TkJ\ 



(104) 



which can be substituted back into ean. (|100l) to eliminate 
A. Defining a = c\-\- 3c2 -I- C3, the resulting zeroth-order 



1 



^-TKa\W ^-TM^a^ = 



SttGoo^ 



E' 



6 3 

- [1 - 2aTK\ U^ -Tk. 1 - \^olT_ 



v^A" = zn 

and the kinetic scalar /C simplifies to 



tC = 'dc 



AP 



(105) 

(106) 
(107) 

(108) 



The first two equations above can be written in the form 
of eqns.® and ([6]) with the identifications 

a^X = iFKofhL^ - -a^TM^ (109) 

a^Y = -2FKaH^ + ^a^J'M'^ - aiFKU) (110) 

For simplicity we will treat only the linear Einstein- 
Aether theory, in which J^{K,) = K,. In this particular 
instance eans. (|105p and (I106P can alternatively be rewrit- 
ten (making use of eqn. (jl08p such that they differ from 
the equivalent expressions in GR only through a modifi- 
cation to Newton's gravitational constant: 



Go 



Go 

1 - ^ 

^ 2 



(111) 



The situation in Einstein- Aether theory is similar to the 
toy model we considered in §IVD| with the instance 
g — a/2. We therefore expect the perturbations of C/^i/ 
to look similar to egns.ip^- PS)) . However, we have now 
introduced an extra field into the theory that was not 
present in the toy model, and this will give rise to new 
terms. Recall that when we applied a gauge transfor- 
mation to eqn. p^ any terms produced by gauge- variant 
quantities cancelled due to the background equation (|40l) . 
From this we can deduce that any new terms introduced 
by the vector field must be explicitly gauge-invariant. 

Now we wish to consider linearly perturbed Einstein- 
Aether theory; we write the perturbations to the vector 
field as 



A^ = ( 1 + Z, -d'Q 
a 



(112) 



Taking the linear perturbation of the constraint A^A^^ = 
— 1 we find that the perturbation to the 0th component 
of the vector field is tied to the metric perturbations by 
Z = S = — '5, the last equality being true in the confor- 
mal Newtonian gauge only. The perturbed equation of 
motion for the vector field is: 



ci(Q + k^Q + 2'HQ + 2'H^Q + * + $ + 2%^ 



+ (3c2 + C3) \k^Q + 2n'Q-- 
' a 



$ + -H* =0(113) 
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One can then find the hnear perturbations of U,,,i, . This 
is a lengthy but straightforward exercise [zl, UM ■ Using 
eqn. (jll3| ) to remove second-order derivatives from Uq^ 
the resuhs are: 



Ua = {ci-a)k"HQ + cik^Q + ^] -3an 



$ 



Ue ^ a 



<^ + n^ 



+ (Cl +C2+ C3)k'^Q 



at 



(q + 2'HQ 



Up 



C/e = -(C1+C3) Q + 2HQ 



3a $ + 2-H$ + (2-H + W)"^ + U^ 



(114) 
(115) 

(116) 
(117) 



It is straightforward to read off the coefficients introduced 
in eqn.(|49|): 



(118) 



ao = Hici ~ a) oci = cik 
Ao = Fi=Q 



and similarly for the remaining coefficients Pi, . . .Ei and 
Bi, ...Ki ofeqns.(|50 )) - ((52| . Note from eqn. (fTT5)) that the 
d.o.f. X (= Q) is not dimensionless, so neither are ag and 
ai. 

The Poisson equation in Einstein- Aether theory is: 

-2fc2$ = 8^TGaa'^^p^Ai 



fc^ •HQ(3ci + 2c3)H-ci(gH-*) (119) ^siip 



Consideration of some special cases of the Ci should help 
us gain some understanding of this expression. Firstly 
we note that the slip between the metric potentials is 
sourced by spatial perturbations to the vector field, so 
in the smooth limit we recover $ = 4'. In this limit the 
only modification to the Poisson equation is through a 
constant rescaling of Newton's constant: 



~2fc2$ 



1 



Gq 

£1 
2 



' J2 p^^' 



(120) 



With Q = all explicit traces of extra fields disappear 
from the components of U^u (though we ought to re- 
member that the perturbation to the time component of 
the vector field still remains, 'disguised' as the metric 
perturbation $). We note immediately that the con- 
stants renormalizing Newton's gravitational constant in 
the background and linear-order equations are not gener- 
ally the same, being a and — ci respectively. If we make 
the choice a = — ci, this special case of Einstein- Aether 
theory can be compared to the toy model considered in 
i jIVDI In that example the background gravitational field 
equations were modified but no extra fields were present, 
and we assumed that the same quantity renormalized G 
at both unperturbed and linearly perturbed order. In- 
deed we find that with the choices g = Bq = Q, a = —ci 



and V = V = (to recover the conformal Newtonian 
gauge) the cqns. (H^ - P5|) reproduce eqns. ([m| - pT7| in 
the limit Q = 0. 

Another case of interest is the choice ci+c^ = 0, C2 = 0. 
This causes the first part of the stress-energy tensor of 
the aether (line 1 in ean. (|100p ) to adopt a form akin 
to the Maxwell tensor of electromagnetism [8(J]. Un- 
der these conditions the gravitational slip again vanishes. 
Note that this differs from scalar-tensor and f{R) gravity, 
where the slip could only be zero if the extra fields (treat- 
ing the scalaron as an extra field) were unperturbed. The 
Poisson equation now becomes: 



k^Cl 



2fc2$ = 8n^a' Y. P^^^ -T^{Q + ^Q 



(121) 



In the quasistatic limit the Q term can be neglected. 
However, the 'electromagnetic condition' has prevented 
us from obtaining a relation between Q, $ and ^. So 
in this special case we are unable to package the modifi- 
cations to the Poisson equation as a modified Newton's 
constant, even in the quasistatic regime. This is an inter- 
esting result, because thus far we have always found this 
to be possible. Indeed, for choices other than ci = — C3 we 
can write the slip relation in the usual format $ — TysUp^, 
with slip parameter: 



1 



, ^2 Cl + C3 



(-f) 



fc2 Cl + C2 + C3 

The Poisson equation is then: 
-2fc2$ = STrGoa^ ^ p,Ai X 

i 

Cl 3ci + 2c3 

2r/siip 4(ci -f C3) 



(122) 



(123) 



1- 



risMp 



As before, the term in square brackets acts like an evolv- 
ing Geff. 

Whilst Einstein-Aether gravity has a similar qua- 
sistatic form to scalar-tensor gravity and f{R) gravity, 
its large-scale behaviour is distinctly different. In the 
limit fc — ?► eqn. (|119p reduces to the standard Poisson 
equation of GR. This is partly due to the absence of a 
term like (5($ -|- H*), such as occurred in the scalar- 
tensor and f{R) cases. Terms like this originate from 
couplings between new degrees of freedom and the cur- 
vature scalar R in the action. A direct coupling of this 
type isn't present in the Einstein-Aether case - instead, 
the coupling between the aether and the metric is en- 
forced through a Lagrange multiplier. However, it is not 
true that Einstein- Aether theory reduces to GR in the 
super-horizon regime, because the slip relation still has a 
non-GR form (eqn. (|117p ). 
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D. DGP gravity 

Dvali-Gabadadze-Porrati (DGP) gravity [5j] has re- 
ceived much attention over the past decade. The model 
considers our (3-|-l)-dimensional spacetime to be a hy- 
persurface (brane) embedded in a five-dimensional bulk, 
with matter fields confined to the brane but gravity free 
to propagate into the extra dimension. There are two 
branches to the theory, arising from a choice of sign ac- 
companying a square root in the Friedmann equation 



H'- 



eH SttG 



E- 



(124) 



where e — ±1. Note that we have returned to using 
physical time rather than conformal time here in order to 
make the late-time behaviour more explicit. The param- 
eter re in this expression defines a crossover scale below 
which gravity is effectively four-dimensional, but above 
which five-dimensional effects become important. It is 
determined by the ratio of the four- and five-dimensional 
gravitational coupling constants: 



background, as specified in eqns.® and ([6]) [l| 






Y 






(127) 



where t denotes physical time. Effective 4D field equa- 
tions are obtained by projecting onto the brane [86| : 



^"^^G^i, + K^g^i, = KjII^^ - E^i, 



(128) 



Here g^u is the induced metric on the brane. The tensor 
n^^ (not to be confused with the scalar pressure pertur- 
bation pU) is given by: 



n„ 



where 



rri rjia 
4 '^" "^ 



(129) 



-^ fif — ^ fiiy 2 ^M^-' 



(130) 



2kI 



(125) 



In a CDM-dominated universe we have H — >■ e/rj, at 
asymptotically late times. The choice e = 1 corresponds 
to a universe that accelerates without the need for a cos- 
mological constant, although the crossover scale must 
still be fine-tuned to fit current supernovae data, with 
Tc ~ Hq^ [55]. Alas, the self-accelerating branch of DGP 
has been effectively ruled out as a viable theory due to a 
ghost-like instability [81-83] , and its failure to fit multiple 
observational data sets simultaneously - a feat achieved 
by ACDM [SJ]. Nevertheless, we propose to describe 
the relations between DGP gravity and the PPF formal- 
ism as a typical model for cosmological perturbations in 
higher-dimensional theories. 

We will adopt the Gaussian Normal longitudinal 
(GNL) gauge, in which the brane remains unperturbed 
at the hypersurface y = and we recover the familiar 
4D conformal Newtonian gauge on the brane only. The 
action for a simple DGP model is: 

1 



S ^ -^ (fX^/^^g (5)i?_2A 



J^^^^R + Ci 



+ / d^x V^ 



ZK^ 



(126) 



where ^^'g and ^^'R are the five-dimensional metric and 
curvature scalar of the bulk; X^ and x^ denote five- 
dimensional and brane coordinates respectively. Cb is 
the Lagrangian of the (brane-confined) matter fields and 
possible brane tension A, which is related to the 5D and 
induced 4D cosmological constants |85| . 

It will be useful to write the modifications to the Fried- 
mann and Raychaudhuri equations in the form of an 'XY' 



E^^ is the projection of the 5D Weyl tensor onto the 
brane: 



E^,^, =(^) CicnnAn-^gjlg:; 



^c B n 



(131) 



By supplementing the Bianchi identities with the Codazzi 
equation and Israel junction conditions [83] one finds that 
the matter energy-momentum tensor is separately con- 
served from li^y and E^v [Sa l: 

V^T„, - 






v^^,,,, = K^v^n^, 



(132) 
(133) 



where V^ is the covariant derivative associated with the 
induced 4D metric on the brane. 

Gomparison of eqns.(Il]) and p28p enables us to 
straightforwardly write down the tensor U/j,^ for DGP 
gravity: 



f/,„. = -A 



4.9m'' 



'^sll/ii^ Efj,u K^l f^i^ 



(134) 



Now the key question is: are we able to perturb this 
expression for C/^i, to obtain UA,Ue,Up and Us in the 
same way that we have done for scalar-tensor gravity, 
f{R) theories and Einstein- Aether theory ( WYMWU^ l 
The answer is a qualified yes and no. We are able to write 
down expressions for these quantities and determine the 
coefficients laid out in ifVl However, the perturbed com- 
ponents of Uf^i, contain four new d.o.f. arising from the 
perturbations to the projected Weyl tensor. Hence, un- 
like the previous three examples presented in this section, 
the system is not closed using just the Bianchi identities. 
We cannot obtain the closure relations needed to solve for 
these new d.o.f. using the purely 4D formalism developed 
in this paper. 
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Let us elaborate. We will label quantities related to 
the Weyl tensor by the letter E, as we shall see it plays 
the role of the effective fluid discussed in fjv] Since the 
Weyl tensor is traceless by construction (and we are as- 
suming isotropy amongst the three spatial dimensions of 
the brane) we deduce that the Weyl 'fluid' must have a 
radiation- like equation of state, we = 1/3. It is common 
practice to neglect the small contribution of the Weyl 
fluid to the cosmological background, and define its com- 
ponents at the perturbed level: 

-E" 



- ^^E 

'E° = WSe 
E] = -fiES) + D)aE 



(135) 



The perturbations of eqn. (|134p are then: 

1 

— / 
6 



TT "42-2 J 2 22r 

ua = Tj'^s'^ p — a liE — ^40 po 



(136) 



Up = -4a^[{p + P)p6 + p^U]-a^PE-3Kla'^pU 



1 
12 



U^ = -^Aty(p + 3P)(p + P)S-aVB 



Kia^p + P)^ 



The quantities marked by tildes are components of T^^ 
defined in ean. (ll30p . and the unsubscripted quantities 
refer to ordinary matter components, as in previous sec- 
tions. By using eqns.(|4]) and (|130p we can rewrite the 
above expressions as: 



the previous examples, we now have a vector of three 

extra d.o.f.: x = {M-BjQ-BjCtb}- 

The two perturbed components of eqn. (|133p are [l[ : 

PE+^npE-'^^QE = (139) 

1 / 7/ 

Qe + ^UQe - T,P^E + (1 + We) Pe + 3 — Ob 
3 \ a 

\E+'-^^^[i2 + 3wE)-^-n]-0 



1 + SwE V 3 



(140) 



Now the difficulty is apparent - we have two equations for 
the three variables p,E , 0_b and ctb . We could eliminate 
fiE or Qe from the above equations, but the anisotropic 
stress of the Weyl fluid remains a free function. If we stay 
entirely within the bounds of a 4D formalism it must 
be treated as an additional source in the 4D effective 
Einstein equations. Some authors have obtained a closed 
system on the brane by setting as to zero, e.g. [88{, but 
this will not be the case is general. 

It is not surprising that ctb appears as a free function in 
the perturbed effective Einstein equations. Gravitational 
waves propagating in the bulk spacetime contribute to 
aE when they impinge upon the brane, and these cannot 
be fully described by brane-bound perturbations. How- 
ever, by tackling the full system of perturbations in the 
bulk it is possible to express the impact of the gravita- 
tional waves on the brane in terms of other brane-bound 
quantities. Mukohyama [831 ^^^ shown that all 5D met- 
ric, matter and Weyl-fluid perturbations can be related 
to a master variable fl which obeys a partial differential 
equation in the bulk. This is the master equation [5Jj: 



Ua 



Ue = 



(Ea + o^PLe) 
iEs+a^OE) 



2rlX 
3 



Up 



2rlX 



Ef 



' 2 
-a Pe 



Xrl 

r^— (1 + We) (Ea -Ua + cl^Pe) 



2Xrl - 3 



C/s = - 



rl{X + 3y) 



{Et. + aVfi) 



(137) 



Well below the crossover scale we recover GR, as C/i — > 
if Tc --?> 00. Using the definitions for the i?i, we can then 
put the above expressions into the form of eqns.(j49 |) - ([52|) . 
For the V a component the non-zero coefficients are: 



^0 



IX 



Ax 






v(l) - 



974 

r^X 

3a^ 

2r2X 



(138) 



where we have dropped the factor of k^ in the first term 
of eqn. p^ in order to keep ap dimensionless. Unlike 



d_ 
di 



1 on 

na^ dt 



d'd. 






nn' 



= 



(141) 



where primes denote derivatives with respect to the bulk 
co-ordinate y, and we have assumed a Minkowski bulk. 
The perturbations of the Weyl fluid are given in terms of 
fiby H: 



Pe 



e. 



aE = 



3a5 

3a* 
1 



-n 



b 

'dt 



Hfl 



sS-SH^-^ + Kn-^'-^n' 



dt^ 



dt 



H dt 



(142) 
(143) 
(144) 



where a subscript b denotes evaluation on the brane and 
k is the 3-momcntum on the homogeneous background. 
Eqns. (|14ip and (|144p arc the expressions we need to close 
the system of perturbations. We have introduced an ex- 
tra variable, il, but this has been compensated for by the 
addition of two new equations. A solution of the master 
equation is then sufficient to fully determine the system 
of perturbations. The fact that all perturbations can be 
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related to a single scalar is a result of the high degree 
of symmetry present in this system (a maximally sym- 
metric brane in a maximally symmetric bulk); fl is not 
generally believed to have a physical interpretation. 

But we are not quite out of the woods yet. Though 
the system of equations is now closed, it is not closed 
on the brane alone. Eqns. (ll41[) and (|144[) depend on 
derivatives of the master variable normal to the brane, 
and hence in general knowledge of the bulk perturba- 
tions will be required to solve the system. Koyama & 
Maartens [901 neatly sidestepped this issue by consider- 
ing perturbations in the small-scale, quasistatic regime 
of a Minkowski bulk. This enabled them to make the 
approximation \Hfl'\ ^ k^ft/a^ so that the troublesome 
terms can be neglected, as well as time derivatives. Un- 
der these assumptions one finds that the energy density 
and anisotropic stress perturbations of the Weyl fluid are 
simply related by ^e = 2k^aE- In the quasistatic regime 
the master equation then has an analytic solution [90]: 



n^C{\ + Hy)-^ 



(145) 



where C is a constant, and regularity of the bulk per- 
turbations has been used to eliminate a second possible 
solution. Solving for the Weyl fluid and metric pertur- 
bations, the modified Poisson equation and slip relation 
are: 



Go 



e{^-^) = S.^a^J^p^A. 



(146) 
(147) 



where /3 = 1 - 2rcH 1 



1 dH 
ilPltt 



It is easy to see that our PPF function g should be iden- 
tified with (1 - 3/3)"\ Combining eons. ([US)) and pl7)) 
to eliminate matter terms, we obtain: 



$-* 



1-3/3 



.$ 



(148) 



From this we can deduce the PPF function 
(■ = 2/(1 — 3/3). Alternatively, in the (Q,r7siip) format 

we obtain: Q= (l - ^) ,?7,iip = (3/3 - l)(3/3 + 1)-^. 

For scales larger than the quasistatic regime the deriva- 
tives normal to the brane cannot be ignored. Sawicki, 
Song and Hu were able to evolve large-scale modes nu- 
merically by implementing a scaling ansatz [9 If : 



n = A{p)aPG{x) 



(149) 



where x is the distance from the brane in units of the 
causal horizon and p is a constant. The authors took 
a trial value for p and solved the system of equations 
iteratively to obtain subsequent corrections. They recov- 
ered the quasistatic solution of Koyama and Maartens 
and were able to calculate the behaviour of horizon-scale 



modes. However, the labour involved in obtaining these 
solutions is decidedly non-trivial. If we wish to constrain 
general classes of gravitational theories simultaneously 
we need a method that does not require detailed numer- 
ical evolution for each theory individually. This is the 
motivation behind the GDM-based approach we put for- 
ward in S|Vl In the case of DGP gravity the new d.o.f. 
are already written in the form of fluid perturbations, 
and we suggest that other theories may be amenable to 
a similar treatment. 

Returning briefly to the quasistatic limit of eans. (|146p 
and (|147l) . we note that the master- variable route to ob- 
taining a modified Poisson equation and a slip relation 
was decidedly different to that taken in WI AIIVl"Cl So 
it is interesting to find that we have reached the same 
conclusion in all four cases: the usual two-function PPF 
parameterization works well in the quasistatic limit, but 
for scales larger than this it no longer captures the full 
behaviour of the theories studied here. 

We should not be particularly surprised by this conclu- 
sion. Given the different physical mechanisms employed 
by the four theories studied here, it is quite a remarkable 
feat that they can be mapped onto a common frame- 
work at all, in any limit - let alone a framework simple 
enough to express the departures from GR using only 
two functions. It is not unexpected to flnd that the cor- 
respondence between parameterization and theory does 
not hold perfectly for all scales. 



VII. DISCUSSION 

The wide array of modified gravity theories now 
present in the literature renders the individual testing of 
theories time-consuming and impractical. There is also 
the possibility that none of the theories currently under 
consideration are correct, or that General Relativity re- 
mains valid for all length scales and environments. The 
Parameterized Post-Friedmann formalism provides a use- 
ful framework for constraining deviations from General 
Relativity without recourse to a specific underlying the- 
ory of modified gravity. 

In this paper we have shown how the constraints of 
energy conservation, gauge invariance and second-order 
field equations can be used to reduce the components of a 
purely metric theory to two free functions. This process 
has highlighted some consistency conditions that arise 
when adopting a parameterized form for modifications 
to GR. We have put forward a general structure for the 
perturbed Einstein equations that should be applicable 
to any single-metric theory for which the gravitational 
action is built from only curvature invariants. However, 
when additional degrees of freedom are introduced into a 
gravitational theory (such as extra fields), our ability to 
make precise statements about the form of the modifica- 
tions to the Einstein equations is reduced. 

Though a rich zoo of underlying physical mechanisms 
has been put forward, ultimately we are concerned with 
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their effects on observables. The observables themselves 
(galaxies, CMB photons) are predominantly controlled 
by the evolution of matter perturbations. There are only 
a limited number of ways that a theory of modified grav- 
ity can affect the matter perturbations: for example, by 
changing the strength of gravitational coupling or the re- 
lationship between the metric potentials $ and ^. But 
these effects are, to a certain extent, degenerate with the 
presence of a second fluid [43 ■ A reduced gravitational 
coupling could be mimicked by pressure support from 
another fluid. A non-zero gravitational slip could be in- 
troduced by a fluid with non-negligible viscosity. Hence 
we suggest that it may be possible to treat the extra d.o.f. 
of a theory by an effective fluid; Generalized Dark Mat- 
ter (GDM) provides one method of doing this in the case 
of spin-0 perturbations. An effective fluid-type approach 
of this kind is already present in DGP. 

Such treatment is not strictly necessary when two or 
fewer new degrees of freedom are present. But if we wish 
to keep our parameterization as general as possible, then 
we need a method to deal with more than two degrees of 
freedom. The GDM approach provides closure relations 
for such cases. The GDM parameters {we, c^g, Cvis} can- 
not be assigned the direct interpretation they possess for 
real fluids. But if constraints on these parameters favour 
bizarre-seeming values that would be unphysical for a real 
fluid than this could indicate their origin to be modified 
gravitational laws and not dark energy. 

To date, the mapping of theories onto the PPF frame- 
work has only been computed explicitly for a small num- 
ber of cases. In this paper we have added Einstein- 
Aether to this collection. In the quasistatic limit the 
two- function parameterization has always been found ap- 
plicable, but it is possible that there are some classes of 
theories which cannot be reduced to such a simple for- 
mat. In addition, we have found that on horizon scales 
the parameterization of the modified Poisson and slip 
relations no longer matches onto the underlying expres- 
sions exactly. These large scales are important for accu- 
rate calculation of how the ISW effect and matter power 
spectrum are affected by modified gravity theories. So 
should we give up on the goal of a unified parameteriza- 
tion? 

Not at all. Consider the analogous problem in dark 
energy; the equation of state is commonly reduced to 
just two numbers via the CPL parameterization [93.l93|: 



w{a) = vjQ + Wa{l — a) 



(150) 



Whilst it is unlikely that a physically-motivated dark 
energy model will map neatly onto this expression, 
ean. (jl50p provides a useful way to obtain constraints on 
the expected behaviour of dark energy. Of course the 
approach is not ideal - if the equation of state were to 
behave in a radically different way to our expectations 
then it would not be adequately described by {wQ^Wa}- 
But the form of ean. (|150p has physical motivation, and 



gives us a way to tackle large classes of models without 
specialising to a particular theory. 

The PPF formalism should be viewed in a similar way. 
Egns. iP^ and pip may not match up exactly to all the- 
ories on all scales, but they do provide a phenomeno- 
logical way to search for the approximate signatures we 
expect modified gravity to leave. Similarly, our proposal 
of mapping additional degrees of freedom onto a GDM 
framework may not hold exactly for all possible theories. 
But whilst we have not yet reached an ideal parameteri- 
zation of modified gravity, we believe that the approach 
outlined in this paper is likely to have a wider range of 
applicability than most of the forms currently in com- 
mon use. If, when combined with the next generation of 
experiments, we fail to make progress on breaking the de- 
generacy between dark energy and modified gravity, then 
we will need to rethink our tools. For the present, use of 
the parameterizations described here is a justifiable sim- 
plification, provided that we are not lulled into thinking 
that their correspondence to underlying theories is exact 
in all situations. 

If we are to use such a phenomenological methodology 
we need to choose a parameterization that closely mod- 
els our expectations. Small differences in the structure 
of the modifications to Einstein's equations can lead to 
significant effects, independent of the ansatz being used 
for the PPF functions. We will demonstrate these differ- 
ences and how they affect the constraints obtained in a 
future work [4d|. 
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Appendix A: Gauge-Invariant Equations for 
Scalar- Tensor Theory 

In this Appendix we display explicitly the gauge- 
invariant form of the modifications to the perturbed Ein- 
stein equations in the case of scalar-tensor gravity. At 
the level of the homogeneous, background universe the 
modifications are contained in the diagonal components 
of a tensor U^^, see eqn.Q. The perturbed components 
of this tensor appearing in eqns.(rT9l)- ([22)) are denoted 
by Ua,Uis,,Up and C/s as given in egns. pH)) . In WIAI 
these components were derived in the conformal Newto- 
nian gauge for simplicity. Using the notation introduced 
in mill the expressions in a general gauge are: 
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C/s = Ss(l- 



# -W^ + ^+^? -2hU(0)44+^^ 



In these expressions the scalar field perturbations Scf), 5(j} 
and 5(j) appear only in gauge- invariant combinations with 
the variables V and S. T is also gauge-invariant. The 
last terms in Ua^Uq and Up are not gauge-invariant, 
and neither are Ea,Eq and Ep. But we find that the 
additional terms produced by these parts under a gauge 
transformation cancel by virtue of the background equa- 
tions (|62p and (|63)) - see Table Hill for the transformation 
properties. Note that all terms in the spatial, traceless 
Einstein equation are individually gauge-invariant. 

So the perturbed Einstein equations with the addi- 
tions above remain fully gauge- invariant, as they must 
do. However, on first inspection it looks like we have 
violated the constraint of having second-order field equa- 
tions. Recall that this constraint restricts Ua and Uq 
to contain at first-order time derivatives at most, due to 
eans. pT]) and P^ . It is easy to see that the (j) terms 
present in Uq cancel, but in Ua this is not explicitly 
obvious. In addition, Ua appears to contain a second- 
order time derivative of the scale factor coming from the 

$ term (see ean. ([T5|) ). To show this is not a problem we 
will need to use the zero-order Bianchi identity, eqn.([7|). 
In the case of scalar-tensor theory this becomes: 



-3 -H 



obtain the alternative form: 
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(A6) 



We see that the terms containing second derivatives of 
(f> and the scale factor have both cancelled. The above 
expression is now demonstrably first-order in time deriva- 
tives, but its ability to yield gauge-invariant Einstein 
equations is no longer obvious. 



Appendix B: Constraint Equations 

In the case of unmodified background equations, con- 
straints are obtained by applying the Bianchi identity 
to a [/-tensor of the general format shown in eans. ([23) )- 
(1^51) . It is possible to deduce the general structure of 
these constraint equations and use this as a shortcut. 
The constraint equations arising from the coefficient of 
$^") can be generated from the following formulae: 



where we have made use of eqn. (|62l) . Using this to sub- 
stitute for the ui(j) term in eqn. ljAll) and simplifying we 



[Bl] An + kAn-i+nAn- 
[B2] Bn + kBn-i+2HB„ 



kB„ 
1 



HCn = 



- Ki^Yl 



TkD„ 



(Bl) 
(B2) 



where Bl indicates ean. pl|) . B2 indicates eqn. p2|) . The 
constraints arising from the term T^"' are analogous: 



[Bl] Fn + kFn^i + HF,, + kin +nJn=0 



(B3) 



[B2] L + kIn-i+2nin-^kJn + ^kKn = (B4) 
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Note that for an Nth-order theory, Fjv_2 and /jv-2 must 
be set to zero to avoid picking up (N-l-l)-th derivatives 
of the scale factor. 
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